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In order to study the dynamic stability of the system, having a precise 
dynamic model including the energy generation units such as generators, 
excitation system and turbine is necessary. The aim of this paper is to design 
a power stabilizer for Mashhad power plant and assess its effects on the 
electromechanical fluctuations. Due to lack of certainty in the system, 


designing an optimized robust controller is crucial. In this paper, the 


establishment of balance between the nominal and robust performance is 
Keyword: done by the weight functions. In the frequencies where the uncertainty is 
Hb/H../u Controller high, in order to achieve to the robust performance of the controller, u 

analysis is more profound, otherwise, in order to achieve to nominal 
LMI i a8 performance, robust stability, noise reduction and decrease of controlling 
Multi-Objective Control signal, the impact of the controller H2/H» is more profound. The results of 
Power System Stabilizer the simulation studies represent the advantages and effectiveness of the 
suggested method. 
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1. INTRODUCTION 

The operation of the power system is usually limited to the boundaries of the dynamic stability, 
which is far from the limits of thermal stability. Nowadays, to improve the attenuation of the oscillations with 
low frequencies, in most cases, power system stabilizers (PSS) are widely used in power systems. Designing 
a proper stabilizer could decrease the limitations due to dynamic issues and help the system to get closer to 
its nominal capacity. These stabilizers are usually designed according to the single machine, infinite bus of 
the system in a definite working point. Therefore, it is possible that the stability of the system to be threaten 
by changes in the parameters or equilibrium points of the system. In this paper, the system stability, in the 
other words the stability of power system’s ability against the change of parameters are checked. As a result, 
three main controlling goals will be obtained: strengthening the closed loop system, lower cost designing 
strategies and improving the transient response. H2 control is used in the transient stage for presenting a fast 
dynamic response to minimize the response energy of the impulse. While H.. control is used in the stable 
stage for reduction in the disturbance and protecting against tracking; which in turn would guarantee the 
robust stability. In order to achieve to optimal performance, taking into account the effect of uncertainty 
during system design period is required. But on the other hand it can lead to severe restrictions on the 
controller that sometimes makes it an infeasible problem. So far, various studies have been conducted on the 
stability and controlling of power systems. 

The first formulation of the H. control problem was performed in 1981 by Zames. To date, large 
numbers of researches have been performed for study of the robust control, the H2 control and He control. 
Doyle has analyzed the state space by using H. and H3 standard form and it’s solving. The conditions of 
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solving problem and its solution using Hamiltonian matrix introduction are the highlights of this paper [1]. 
Also Doyle as well as a tutorial overview in the fractional linear transformations (LFTs) and the value of the 
unique structure, p, and linear matrix inequalities (LMIs) in the solution of LFT problems has offered [2]. 
H/H», were combined by Rotea in this way, two important approaches were suggested, 1) optimal control 
limit of H2 and Hə (actually constrained optimization), and 2) at the same time optimal control of H2/H. [3]. 
Lanzon in his PHD thesis chooses the weight functions in u and He design [4]. Many of the power stabilizers 
proposed for systems of the single machine are not able to resolve the interaction problems; while some of 
the multi variable stabilizers are also lacking suitable robust stability. Studies on the stability are mostly 
conducted on two transient and steady states. At operation condition, a power system is in its permanent state 
[5]. When performance is in the permanent state, if a sudden change happens, the system will go toward 
the disturbance. 

Investigation of the classic stability [5], the optimization method with the help of pareto multi- 
objective [6], the method of adaptive control [7], the nonlinear controller [8], using the parameters estimation 
[9], robust controller H2/H. [10], the pole placement and application of the linear matrix inequality [11], 
fuzzy and Neural network control [12], and Evolutionary algorithm [13] are among the works which had 
been done. The problem of closed-loop identification of the Heffron-Phillips model parameters is of practical 
importance since the data used for identification can be gathered when the machine is normally connected to 
the power system [14]. In this paper, at first the power system was modeled. Then, the problem was 
introduced. In this paper, at first, H2/H. controller was investigated with a new insight along with the new 
controller, u; and then these two different controllers were combined via the use of the weight matrices. 
Solving this problem would be possible by application of the linear matrix inequality. The results indicate 
that, the goals of H»/H./u combination, including elimination of the perturbation effect, reducing the 
controlling signal and accounting for the uncertainty during the system’s functionality investigation, were 
properly realized. 


2. POWER SYSTEM MODELLING 

The stabilizers of the power systems are designed with the aim of improving the attenuation of the 
low frequency oscillations of the system, based on the single machine, infinite bus model. The power system 
stabilizer is a traditional and economic controller whose aim is to increase the dynamic stability of the power 
system. By creating the damping electrical torque, the stabilizer of the power system will improve the 
deviations of the rotors rotations. The mentioned equipment also optimizes and tunes the exciting voltage, by 
creating the suitable voltage. The power plant of Mashhad city is located at the eastern part of the city at the 
beginning of Sarakhs Boulevard. This is the oldest power plant of Khorasan province and has 8 electricity 
generating units, 4 of them are working with steam and the other 4 ones are gaseous. The steam units consist 
of two ELIN and two SKODA units, and the gaseous units include two BBC units and two ALSTOM units. 
This power plant was established in 1964 and started its work in 1968. The exciting system of ALSTOM 
gaseous units of the power plant of Mashhad are classified as the static type. Feeding of such exciting system 
is done via power voltage transformer and three current transformers [15] with the capability of being 
saturated. The controller part of the stimulation system of ALSTOM gaseous units includes 3 main control 
modules. By elimination of the three controlling modules, in order to attenuate the oscillations, the power 
system stabilizers could be applied. In studying the dynamic stability of the power networks, and also in the 
cases where the changes and disturbances of the network are mainly partial and slow, the linear generator 
model could be employed. In order to consider a synchronous generator, we use 3 rd order synchronous 
generator model called Heffron-Philips model [11]. This model contains 3 state variables: Aw,, Ad, AEg. 

Considering the exciter model will lead to the introduction of the fourth state variable AE». In this 
model, governing differential equations are linear around operating point. Figure 1 shows block-diagram of 
linear mode of Heffron-Phillips model along with exciter and AVR Regarding the generator parameters 
Hefron-Philips coefficients could be obtained by (1) [16]: 


' 


E,E_ cosd X,-X! X,+X, Xqy-X 
Kica + Epsinð „Ko == i Ka = —*E, sind 
Xg +X Uae Xa; X, +Xı (1) 
KEX, E,X, cosd x’ XE 
ee eee PEA dE sind, Kg =— 
me a (ARV, Oa, (Xe +X )V, 


State space of equation of Figure 1 shows in (2). 
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Figure 1. Heffron-Phillips model 


3. PROBLEM STATEMENT 
3.1. H2/H. Controller 

Existence of uncertainty created due to an uncertain and erratic input (noise and disturbance) and 
Un-modeled dynamic cannot be described completely and precisely as a true system by a mathematical 
modeling. On the other hand, a true system should contain the following important objects: robust stability, 
robust and nominal performance, settling time, maximum over shoot and etc which try to gain these 
objectives of the controlling problem [4]. The type of uncertainty is another important factor in the system 
analysis. Consider additive uncertainty shown in Figure 3. 


Figure 2. M — A Model Figure 3. Additive uncertainty 


=- 
Objective 1: if A = 0 then |-s| < 1 (nominal performance). S = (J + GK) (S is sensitivity function). 
Objective 2: if A # 0 then system is robust stability. M = (1 + KG) 'K, 

if EAG) < jo) = |m] <1 3) 


Objective 3: n is white noise with one PSD (power spectral density). H2 norm, caused due to decrease in the 


controlling signal. Jere, oe (To minimize U; variance with noise input). F(s), R(s) and y(s) are weighting 


function) from Parseval equation and objective 3. Then we have three tasks for controller design 
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FS(K,G) 
yM(K,G) 
RT au, (K,G) 


(FSI, <1, (4) 


| <1), such that, 
a 


(S)M| <1> Thu, 


<1 


o] 


Problem (4) shown in Figure 4. Rotea and Doyle offer two other methods for solve this problem. 
[1]-[3]. A large class of system with uncertainty can be treated as LFT (Linear fractional Transformation). 
LFT model is shown in Figure 3. W: the disturbance signals to the system which won’t be a function of states 
of system, Z: the variable that will be controlled, P: the nominal open loop system, Y: the system measurable 
output. To transform the changed diagram of Figure 4 to the LFT model, we will write the problem in the 
standard form, and then solve it by using of Riccati equation [17]. The (4) LFT model is practicable in Figure 
5 and can be used to design a controller. State space of Figure 5 is written in (5). 


E A 0 00 1x 0 B B T 
ty _ -B,C A; 0 0 Xf if Br -B,D -B,D ee 
i 0 0 A 0 lx jo 0 B 
. y 7 4 Ww 
iy 0 0 0 Apr] [0 BR Bp 
x = Ax +B,W +Bou ge 3 B, | 
zZ=C,x+D,,;,W+Dy,5u = 5 
1 11 12 4 -D;C C; 0 O]fx D; -D,D -D,D T (6) 
y=C2x+DaW+D2u jalo o c olx do o Bl ae 
% 0 0 0 Cell} }0 Dg Dp || 
Ly = 0 0 0 |x 1 -D -D 
zı] [Fs i 7 
Z=| Z, |=| yM a CL 
2 
Z RT i 
3 nU} 


Determining three weight functions, specified in Figure.4, contain special importance. Using robust 
optimal state feedback method for (4) equations. 


Figure 4. LFT Model Figure 5. Graphical model of problem (4) 


3.2. u Controller 
Here we try to assess robust performance of this closed-loop system by using u-analysis. Robust 
performance condition is equivalent to the following structured singular value u test [2]. 


MA <7" vil <r e my vw (6) 


1 
min {(A)|det(J - MA) = 0 


The complex structured singular value siç) is defined as HAc) = } Lower and 


Upper bond of u can be shown to be P(UM) < u, (M) < min &(DMD ') . 


3.2.1. D-K iteration 
Unfortunately, it is not known how to obtain a controller’s achieving path directly to the structured 
singular value test. But we can obtain the lower and upper bounds of u. This method taken here is the so- 
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called D-K iteration procedure. The D-K iteration involves a sequence of minimizations over either K or D 
while holding the other fixed, until a satisfactory controller is constructed. First, for D = I fixed, the 
controller K is synthesized using the well-known state-space H», optimization method. LFT form of Figure 3 
is written in equations (7) [17], [18]. 


P q 0 00 IP 
x=Ax+[0 0 BIW] , |z Ckx+/-I I -D| W (7) 
U y| |-c -I I -DIU 


3.3. New approach: H2/H., p combination 

Now, we tend to synthesize two collectors according to Figure 6. As mentioned before, the 
availability of robust performance causes extreme limitation on the controller, which sometimes prevents it 
from reaching a possible condition. Also, availability of nominal performance means considering operation 
without uncertainty, and it is usual that the essence of uncertainty has decisive effect on the operation. So, we 
tend to balance between robust and nominal performance. W, and W2 are weight functions. Having this data, 
we can determine which frequencies have more uncertainty effect, with regard to the controller effect of p. 
Of course, it is of importance to mention that robust performance contains nominal performance, so, 
controller coefficient of u should be smaller than H/H. controller coefficient. 
Problem 1: Determine W; and W,, in a way that an additive uncertainty system contains robust stability. 


M =(W,K,G+W,K,G+1) (WK, +W,K,) 


8 
[m], <1 i 


PER Ree NEARE E E eee 


— W H, / Hinfinity 
1 Controller + 


|——> —» 


Input wW MU Synthesis| + # Output | 
' 2 Controller 
Controller 


Figure 6. Controller H2/H../u 


3.3.1. Robust optimal state feedback with H2/H., p combination 

We now attempt to follow the analysis of the conditioning of the pole placement problem. 
Researchers shown a number of robust performance indices have been considered in optimization approaches 
for control system design [18]. In robust control using H» optimization, the objectives are expressed in terms 
of the H..-norm of transfer functions. One of the objectives is the following: 


masip ol s| , Where S =[J -(jol —A) | BK a . In this paper we assume that the state of the generalized 
K 

plant G is available for feedback. To be more precise let a state-space description of P (figure 3) is given by 

(LFT Model): 


x= AX +BU+ByW 

Y=X (9) 
U = KX 

ż=(A+ BWK +BW2K3)X + ByW 


The signal W denotes disturbance. The signals U and Y denote the control input and the measured 
output, respectively. Next to gaining Kı by H2/H. and Kz by u analysis, we tend to determine weight 
functions, using linear matrix inequality. 
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Lemmal: (bounded-real lemma) given a constant y > 0, for system, M(s) = (A, B, C) the following two 
statements are equivalent, 1) this system is stable | M(s) || <7 > 2) there exists a symmetric positive 


definite matrix Q, such that: [19] 


Pr T 
A Q+QA QB y C, 
T -1 T 
B, Q -y I D <0 
P y 2 (10) 
C, D; -y I 
Q>0 


Lemmaz2: Consider the feedback system of Figure.3, where G is given by (9). Then, a given controller K is 
admissible and close loop system is robust stability and desired performance if and only if there exists 
W, and W, solving the following LMI problem: 


AQ + BY, +BY, +Y} B+YZB+QA" By act 
By -$I D” |<0 
ca D -pI 


9>0 , Y>0 , %>0 


Where, W, = YOK, W, = YO'K". K, and K, Design with equations 5 and 7 and controller achieves 


Lemma 2, it helps to solve of problem 1. Mashayekhifard et al. presented Robust multi-objective static output 
feedback with H2/H., u combination [20]. 


4. METHODOLOGY 
a. To design the H2/H. for the process with uncertainty. (It helps to select the weighting function properly). 


FS(K,G) 


For H2 /Hə» design can use Rotea and Doyle method. ([3], [8]) or use mM (K,G) < l and obtained K,. 
RT(K,G) 


For F,y and R we use inverse sensitivity function. Or use Automatic Weight Selection Algorithm [4], 
[21]. 
b. To design the u controller for the process with uncertainty (if the process is unstable, at first must be 
stabilize). D-K iteration method can be used to improve the performance of the controller design for the 
system. Peak value of the u (D-K iteration) bound should be less than one, and obtained K, . 
Order reduction method can be used to reduce the order of the K,,K,. 
W,,W, are given with LMI (12) then the robust stability of the system has to be established. 
H infinity norm of W, must be smaller than w, . 
K=W,K,+W,K,. 
This controller (K) has robust stability and desired performance. 


moe o 


5. RESULTS OF SIMULTION 
First H2/H.. controller and then u is designed. After that the order of I+GK was reduced by the help 
of the residual method. Regarding the practical considerations and by application of the inverse of the 
s+2 p-40b+D] 


sensitivity functions, the weight functions were selected with the form of R= l, 
s +200 0.5s +4 
Ti = I . Ki and K: are determined according to the equations 5 and 7, while w; and w2 were defined 
S+ 


regarding equation 12. According to Figure 6 and section 3.A and 3.B, k was designed. The simulations were 
done by MATLAB software and toolboxes of LMI [22], Robust multiobjective control toolbox [23] and u 
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[24] were employed. In the designing process, we used Heffron-Phillips model which is a reduced order 
model. In order to estimate our designs through simulations, we use complete model of power system 
containing synchronous generator, exciter system, governor, turbine, 3-phase transformer, transmission line, 
load and infinite bus. For comparison purposes, we compare the variations of before and after 3-phase fault 
occurring in the middle of transmission line. Three-phase fault occurs at 0.5 sec. and is gone within 0.55 sec. 
In addition, the comparison of the singular values for controlling signals related to three types of design is 
depicted in Figure 7. The results show that the largest amount of the control signal is related to u controller 
and the lowest amount was associated with H2/H... Step response of the closed loop system for the three 
controllers shown in Figure 8. 

Figure 8 indicates that the best function of the controller is for u while H2/H.. shows the weakest 
performance. It could also be noted that since the system has multiple inputs and outputs, the sensitivity and 
weight functions have the matrix form. The results verify the success of combining the robust and nominal 
performance with each other. Reaching to the mentioned objectives with the minimum controlling signal is 
one of the advantages of H2/H../u controller. Most of the robust controllers have high orders and controlling 
signals. But this new approach did well in this regard. H2/H. controller has the order of 7, and pu controller’s 
order in 10, due to use of order reduction method, the order of the H2/H../u controller is 5. For further 
investigation of three controller, the form of the waves related to rotor angle and speed are shown according 
to 2% p.u increase in the input voltage of the system in Figure 9 and 10, respectively. The mentioned figures 
indicate for H2/H../u controller the attenuation rate of 3 s and low oscillation. Variations of rotor speed before 
and after 3-phase fault shown in Figure 11. 


Controlling Signal 


H2/Hinf/MU 
ms AL ma T H2Hin f 
3 MU 


Singular Values (dB) 


Frecuency (rad/s) 


Figure 7. Singular value for controlling signal (H2/H», u, H2/H./u) 


step response of close loop system 


H2/Hinf/MU 
H2/Hinf 


A 
/ Ja 
N 
17 MU 


Amplitude 


Time (seconds) 


Figure 8. Step response of close loop system (H2/H.., u, H/H/u) 
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Figure 11. Variations of rotor speed before and after 3-phase fault 


6. CONCLUSION 

Providing the spare parts and resolving the errors in the excitation system are among the most 
important problems of the old power plants. For this reason, replacement of the control section of the 
excitation system seems necessary. For attenuating the oscillations by controlling the excitation process, the 
stabilizers of the power systems are used. The aim of this paper is to design a robust stabilizer of the power 
system for the power plant of Mashhad city. First the parameters of Hefron Philips model was derived and 
obtained, since there is no certain model of the system in hand, the robust performance is considered. By 
robust performance, it means by consideration of the uncertainty the errors of the system be minimized. In 
order to investigate the robust performance, u analysis was used. Generally, existence of the robust 
performance results in the severe limitations on the controller which is sometimes making it an unfeasible 
issue, and if it could be feasible the order of the controller would become higher and the resulted control 
signal would be increased which would lead to saturation of the actuator. In order to decrease the controlling 
signals, it is needed to use to controllers of u and H2/H. for the performance of robust and its stability. 
Designing the filters or in the other words weight functions have also crucial role in determination of the 
closed loop response. In this content, first, three weight functions were designed for H2/H. controller and 
then two weight functions by LMI were designed for balancing between H2/H.. and u. Due to multi variable 
system of the weight functions, they were plotted in the form of matrix and the singular values. The results 
show that the closed loop was stabilized despite of the existence of uncertainty and has the desirable 
performance. Moreover, the response of the closed loop and controlling signal of the combined controller 
(H2/H.,./u), is between the two other controllers. The angle and speed of rotor verifies the effectiveness and 
advantages of the suggested method. 
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7. APPENDIX 
A. Nomenclature 
Xa direct axis reactance of synchronous machine (p.u) H inertia constant 


X’q direct axis transient reactance of synchronous E» exciter Output Voltage 
machine (p.u) E, voltage proportional to direct axis Flux 
Xq quadrature axis reactance of synchronous machine linkages 
X’, quadrature axis transient reactance of synchronous 6 (t) rotor angle 
machine ©r (t) speed of the rotor 
X. transmission line reactance Tm mechanical/electrical torque 
T'ao direct-axis transient open circuit time constant Iq generator stator current 
Ki to Ke Heffron-Phillips model coefficient V: Teminal voltage of synchronous machine(p.u) 
Ka DC gain of the AVR A Denotes small perturbation in the variable from 
Ta time constant of the AVR steady state value 
Kp PSS gain fy Synchronous Generator 
B. Machine data 
Xa X4 Xi X4 Xe Ta Kı K2 
2.013 0.3 1.76 0.65 0.68 0.53 0.55 1.2 
K; Ky K; Ks Ka Ta Kp H 
0.66 0.7 0.095 0.815 50 0.5 7.1 3.5 
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